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Doubly charmed mesons are studied within a quark model framework. We solve the 
four-quark Schordinger equation by means of a variational approach by using different 
interacting potentials. Our results point to the existence of a stable isoscalar doubly 
charmed four-quark meson with quantum numbers J p = 1+. 
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The potentiality of the quark model for hadron physics in the low-energy regime 
became first manifest when it was used to classify the known hadron states. De- 
scribing hadrons as qq or qqq configurations, their quantum numbers were correctly 
explained. This assignment was based on the comment by Gell-ManrPJ introducing 
the notion of quark: It is assuming that the lowest baryon configuration (qqq) gives 
just the representations 1, 8 and 10, that have been observed, while the lowest meson 
configuration (qq) similarly gives just 1 and 8'. Since then, it has been assumed that 
these are the only two configurations involved in the description of physical hadrons. 
However, color confinement is also compatible with other multiquark structures, like 
the tetraquark qqqq first introduced by Jaffe^ 1 . 

The idea of unconventional quark structures is quite old and despite decades of 
progress, no exotic meson has been conclusively identified. A conclusive evidence of 
its existence seems rather difficult to be obtained with states that could be described 
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by means of the lowest order configurations suggested by Gell-Mann, like the scalar 
mesons or uncharged charmonium resonances. In such cases, the experimental data 
will hardly exclude a theoretical modeM The study of exotic structures that cannot 
be accommodated in the scheme proposed by Gell-Mann should be of primary 
experimental and theoretical interest, because their existence would make manifest 
the contribution of four-quark states to hadron spectroscopy^. Among these exotic 
structures there are two types that could be experimentally accessible in present 
facilities: charged charmonium states^ and doubly heavy mesons. 

Doubly heavy mesons is one of the first scenarios where the existence of bound 
multiquarks was proposed, systems composed of two light quarks and two heavy 
antiquarks (nnQQ). Although they may be experimentally difficult to produce and 
also to detectP it has been argued that for sufficiently large heavy quark mass they 
should be bouncP. Their stability relies on the heavy quark mass. The heavier the 
quark the more effective the short-range Coulomb attraction to generate binding, 
in such a way that it could play a decisive role to bind the system. Moreover the 
QQ pair brings a small kinetic energy into the system contributing to stabilize it. 

In this talk we revise the possible existence of stable doubly heavy mesons within 
the nonrelativistic quark model. For this purpose we solve the four-body Schrodinger 
equation based on a variational method in terms of generalized gaussiani^. The four- 
quark wave function is written as a sum of outer products of color, isospin, spin and 
configuration terms 

\(f>ciSR >= IColor > (Isospin > [|Spin > (g>\R >} JM , (1) 

such that the four-quark state is a color singlet with well defined parity, isospin 
and total angular momentum. Coupling the color states of two quarks (antiquarks) 
can yield two possible representations, the symmetric 6-dimensional, 6 (6), and the 
antisymmetric 3-dimensional, 3 (3). Coupling the color states of the quark pair 
with that of the antiquark pair must yield a color singlet. Thus, there are only two 
possible color states for a QQqq system, 

|Color >= {|3 12 3 34 >, |6 12 6 34 >} . (2) 

These states have well defined symmetry under permutations. The spin states with 
such symmetry can be obtained in the following way, 

|Spin >= \((si,s 2 )Si2, (s 3 ,s 4 )5 , 34)S' >= |(Si 2 S 3 4)S > ■ (3) 

The same holds for the isospin, |Isospin >= (2 3 ,i 4 )/ 3 4 >, which applies only to the 
n-quarks, thus I = 

To describe the spatial part of the wave function we choose for convenience the 
Jacobi coordinates, 

xi = fi - f 2 , 
x 2 =r 3 -r 4 , 

mir*l+ 7712^2 771 3 f 3 +TO4f4 

X3 = , (4) 

mi + m 2 m 3 + mi 
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Using these vectors, it is easy to obtain basis functions that have well defined 
symmetry under permutations of the pairs (12) and (34). The radial wave function 
is expanded as a linear combination of generalized gaussians 

n n 

\R >= J2 ^ ] R\xi,X2i x 3 ) = P (l) R l (5) 

i=l i=l 

where n is the number of gaussians used for each color-spin-flavor component and 
R l depends on six variational parameters, a 1 , b % , c 1 , cf , e l , and / z , one for each scalar 
quantity: x% , x| , Xi ■ a?2, x\ ■ £3 and X2 ■ X3. Once the spin, color and flavor 
parts are integrated out the coefficients of the radial wave function are obtained by 
solving the system of linear equations 

^2^ i) [(R^\H\R^)-E(R ij) \R^)]=0 Vj, (6) 

i 

where the eigenvalues are obtained by a minimization procedure. 

Assuming nonrelativistic quantum mechanics, the hamiltonian will be given by 

i=l ^ l ' i<j=l 

where the potential V(fij) corresponds to an arbitrary two-body interaction. We 
have analyzed the possible existence of four-quark bound states by using two stan- 
dard quark-quark interactions, a Bhaduri-like p ote ntial (BCN)2J and a constituent 
quark model considering boson exchanges Both interactions give a reason- 

able description of the meson and the baryon spectroscopy, a thoughtful requirement 
considering that in the tetraquarks qq and qq interactions will contribute. 

We have performed an exhaustive analysis of the QQqq spectra by means of the 
quark models described above. We have considered all isoscalar and isovector states 
with total orbital angular momentum L < JH[ i n Table [TJ we show the results for 
some selected quantum numbers obtained with the CQC model. Our study draws a 
first general conclusion: once the four-body problem is properly solved the number 
of bound states is rather small if not zero. For the ccnn only the J p = = 0) 
is bound independently of the interacting potential, already proposed more than 
twenty years ago in Ref. [7] 

We have also studied all ground states of the bottom sector using the same 

interacting potentials. It was pointed out in the early 80's that a QQQ'Q' four-quark 

state should be stable against dissociation into QQ' + QQ' if the ratio rriQ/mQ' is 

large enouglff^l. We observed that all bound states become deeper than in the charm 

sector and a new state, J p = + (I = 0) appears, strengthening the conclusion that 

the larger the ratio of the quark masses the larger the binding energy. 

In connection with the interacting potential used, it has been re cent ly analyzed 
- - • • ji 01 

the stability of QQnn and QQnn systems in a simple string model considering 

only a multiquark confining interaction given by the minimum of a flip-flop or a 
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butterfly potential in an attempt to discern whether confining interactions not fac- 
torizable as two-body potentials would influence the stability of four-quark states. 
The ground state of systems made of two quarks and two antiquarks of equal masses 
was found to be below the dissociation threshold. While for the cryptoexotic QQnn 
the binding decreases when increasing the mass ratio mQ/m n , for the flavor exotic 
QQnn the effect of mass symmetry breaking is opposite. 



Table 1. Four-quark state properties for selected quantum numbers. All states 
have positive parity and total orbital angular momentum L = 0. Energies are 
given in MeV and distances in fm. See text for details. 



{St, I) 


(0,1) 


(Id) 


(1,0) 


(1,0) 


(0,0) 


Flavor 


cenn 


cenn 


cenn 


bbnn 


bbnn 


Energy (E 4q ) 


3877 


3952 


3861 


10395 


10948 


Threshold (T = Mi + M 2 ) 


DD |s 


DD* |s 


DD* |s 


BB* |s 


B\B \ P 


A E = E 4q -T 


+5 


+ 15 


-76 


-217 


-153 


Pmm 


1.000 








0.254 


Pmm* 




1.000 


0.505 


0.531 




Pm*m* 


0.000 


0.000 


0.495 


0.469 


0.746 




60.988 


13.804 


0.787 


0.684 


0.740 




60.988 


13.687 


0.590 


0.336 


0.542 


(4) 1/2 


0.433 


0.617 


0.515 


0.503 


0.763 


RMS 4q 


30.492 


6.856 


0.363 


0.217 


0.330 







Let us finally discuss the role played by hidden-color configurations, color singlets 
built by nonsinglet constituents. Besides Eq. ©, there are other two different ways 
for coupling two quarks and two antiquarks in a colorless state, 

[(?l?3)(t?294)] = {|ll3l24). (813824) } 

[(«i«4)(«a« 3 )] - {II14I23), |8i 4 8 23 )} . (8) 

These two basis contain both singlet-singlet (physical) and octect-octect (hidden- 
color) components. It is possible to prove that there is a minimum value for 
the octect-octect component probability of the wave function either in the 
[(Qi«3)(<22^4)] or the [(Qin 4 )(Q 2 n 3 )} couplings: Pgg 3 ' 24 , Pgg ,23 G [1/3,2/3]. Does 
this imply an important hidden-color component in all QQnn, states? The answer 
is no^, because one can express any QQnn state in terms of the singlet-singlet 
component of the [(<3iH 3 )(Q2fi4)] and [(Qifi4,){Q 2^3)] basis. 

This discussion can be made more quantitative. Let us assume that {P, Q} 
and {P, Q} are the projectors associated to two orthonormal basis that are not 
orthogonal to each other, i.e., PP 4>) ^ and PQ \ 4>) 7^ for an arbitrary 
state I 4>) . This would be the case of the two orthonormal basis in Eq. ([SJ . For any 
arbitrary state | = P | + Q \ \&) , the probability associated to P or P will 
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<J | PQ | *) + (* | QP | * 



be given b>EI, 

J "* , «»l' = 2(l-cos=«) 

P '* ,( [«'l» = 2(1 -L»a) K* ^ 1 *) + (* 1 QP 1 *)] ' (9) 

where P = | u) (u \ and P = | it') (u' | and cosa = (it' | it). For a molecular state 
either ^'^([it]) or T 5 '*^^']) would be close to zero while for a compact state both 
will be different from zero. 

We show in Table [1] the meson-meson probabilities for some selected four-quark 
states according to Ref. [14] Unbound states converge to two isolated mesons, the 
lowest threshold of the system, its root mean square radius (RMS) being very large. 
In contrast, bound states have a radius smaller than the threshold and they present 
probabilities different from zero for several physical states, the lowest two-meson 
threshold being contained in the physical four-quark system. Such states would be 
called compact in our notation. When the binding energy approaches the threshold, 
the probability of a single physical channel converges to one, what we defined as 
a molecular state. The cenn J P (I) = l + (0) is therefore a clear example of a real 
tetraquark, the r.m.s. of all Jacobi coordinates being small and with a significant 
probability of different vectors of the Hilbert space, pseudoscalar- vector (MM*) 
and vector-vector (M*M*) components. 
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